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Abstract 

In this paper we calculate the branching ratios (BRs) of the 32 charmless hadronic B c — > AP 
decays (A = a x (1260), &i(1235), K x (1270), #1(1400), /i(1285), /i(1420), hi(1170), hi (1380)) by 
employing the perturbative QCD(pQCD) factorization approach. These considered decay chan- 
nels can only occur via annihilation type diagrams in the standard model. From the numerical 
calculations and phenomenological analysis, we found the following results: (a) the pQCD pre- 
dictions for the BRs of the considered B c decays are in the range of 10" 6 to 10" 8 , while the 
CP-violating asymmetries are absent because only one type tree operator is involved here; (b) 
the BRs of AS = processes are generally much larger than those of AS" = 1 ones due to the 
large CKM factor of IV^/^usI 2 ~ 19; (c) since the behavior for l P\ meson is much different 
from that of 3 Pi meson, the BRs of B c — > A^P^P decays are generally larger than that of 
B c -> A( 3 Pi)P decays; (d) the pQCD predictions for the BRs of B c (Ki(1270),#i(1400))r?W 
and (i^i(1270), -fCi(1400))i^ decays are rather sensitive to the value of the mixing angle 9k- 
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I. INTRODUCTION 



Unlike the ordinary light B q (q = u, d, s) mesons, the B c meson is the only heavy meson 
consisting of two heavy quarks b and c and plays a special role in the precision test of 
the standard model(SM) Moreover, a large number of B c meson events will be col- 
lected with the running of Large Hadron Collider(LHC) experiments and this will provide 
great opportunities for both theorists and experimentalists to study the perturbative and 
nonperturbative QCD dynamics, final state interactions, etc. 

In two recent works [2|, y| , the pure annihilation B c — > PP, PV/VP, VV decays (here P 
and V stand for the light pseudoscalar and vector mesons) have been studied by employing 
the SU(3) flavor symmetry and the pQCD factorization approach 0-0], respectively. 

In the present work, we will study the two body charmless hadronic B c — > AP decays 
(here A denotes the light axial- vector mesons), which can only occur via annihilation type 
diagrams in the SM. First of all, the size of annihilation contributions is an important issue 
in the B meson physics, and has been studied extensively, for example, in Refs. [H, 0,0- 10 



Secondly, the internal structure of the axial- vector mesons has been one of the hot topics 
in recent years 

HM- Although many efforts on both theoretical and experimental sides 
have been made 14j-|20l| to explore it through the studies for the relevant decay rates, the 
CP- violating asymmetries, polarization fractions and the form factors, etc., we currently 
still know little about the nature of the axial-vector mesons. 

In the quark model, there are two different types of light axial vector mesons: 3 Pi 
and l P\i which carry the quantum numbers J PC = 1 ++ and 1 H , respectively. The 1 ++ 
nonet consists of ai(1260), /i(1285), /i(1420) and K 1A , while the 1+" nonet has &i(1235), 
/ii(1170), /ii(1380) and K\b 1 - In the SU(3) limit, these mesons can not mix with each 
other. Because the s quark is heavier than u, d quarks, the meson l£i(1270) and -Ki(1400) 
are not purely l 3 Pi or 1 1 P 1 state, but a mixture of K\a ( 3 Pl state) and Ki B ( 1 Pi state). 
Analogous to rj — rf system, the flavor-singlet and flavor-octet axial- vector mesons can also 
mix with each other. It is worth mentioning that the mixing angles can be determined by 
the relevant data, but unfortunately, there is no enough data now for these mesons which 
leaves the mixing angles basically free parameters. 

In this paper, we will calculate the branching ratios of the 32 non-leptonic charmless 
B c — > AP decays by employing the low energy effective Hamiltonian 2l| and the pQCD 
factorization approach based on the framework of kx factorization theorem. By keeping 
the transverse momentum of the quarks, the pQCD approach is free of endpoint singu- 
larity and the Sudakov formalism makes it more self-consistent. In the pQCD approach 
one can do the quantitative calculations of the annihilation type diagrams directly, which 
can be seen, for instance, in Refs. j3|-0, 0, 0] . 

The pure annihilation B c — > PP, PV/VP, VV decays and B c — > AP decays considered 
in Refs. [3, 0] and in this paper generally have very small branching ratios: at the order 
of 10 -6 to 10 -9 . According to the discussions as given in Ref. [2], the charmless hadronic 
B c decays with decay rates at the level of 10 -6 could be measured at LHC experiments 
with the accuracy required for the phenomenological analysis, while it may be difficult to 



For the sake of simplicity, we will adopt the forms a± and b\ to denote the non-strange axial-vector 
mesons ai(1260) and &i(1235), respectively, in the following section. We will also use K\ to denote 
-£Ti(1270) and i^i(1400) for convenience unless otherwise stated. 
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measure those B c decays if their branching ratios are much less than 10~ 6 . 

The paper is organized as follows. In Sec. [TTl we present the formalism of the consid- 
ered B c meson decays. Then we perform the analytic calculations for considered decay 
channels by using the pQCD approach in Sec. IIHI The numerical results and phenomeno- 
logical analysis are given in Sec. IIVI Finally, Sec. |V] contains a short summary and some 
discussions. 



II. FORMALISM 

In the pQCD approach, the decay amplitude of the two body decay B c — > M\M 2 
(Mi, M2 stand for the two final state mesons) can be written conceptually as the convo- 
lution, 

A(B C -> M X M 2 ) ~ Jd 4 hd%d% Tr [Cit^B^k^M^k^M^Hik^^k^t)] , (1) 

where fcj's are momenta of light quarks included in each mesons, and Tr denotes the trace 
over Dirac and color indices. C (t) is the Wilson coefficient which results from the radiative 
corrections at short distance. In the above convolution, C(t) includes the harder dynamics 
at larger scale than ms c scale and describes the evolution of local 4- Fermi operators 
from raw (the W boson mass) down to t ~ 0(^/Am Bc ) scale, where A = m# c — m&. 
The function H(ki, k 2 , k 3 , t) describes the four quark operator and the spectator quark 
connected by a hard gluon whose q 2 is in the order of ArriB c , and includes the 0{\/K.m Bc ) 
hard dynamics. Therefore, this hard part H can be perturbatively calculated. The 
function $ M is the wave function which describes hadronization of the quark and anti- 
quark to the meson M. In the present work, since the B c meson is composed of two heavy 
quarks b and c, we will take the nonrelativistic approximation form 8(x — m c /mB c ) [22| 
for the distribution amplitude </>_b c (x). For light meson A and P, we adopt the light- 
cone distribution amplitudes directly, which will be displayed in Appendix |A] While the 
function H depends on the processes considered, the wave function $m is independent 
of the specific processes. Using the wave functions determined from other well measured 
processes, one can make quantitative predictions here. 

Since the b quark is rather heavy, we work in the frame with the B c meson at rest, 
i.e., with the B c meson momentum P 1 = (m_B c /v2)(l, 1, Ot) in the light-cone coordinates. 
For the charmless hadronic B c — > AP decays, we assume that the A (P) meson moves in 
the plus(minus) z direction carrying the momentum P 2 {Pz)i an d with the polarization 
vector €2 for the A meson. Then the two final state meson momenta can be written as 

^2 = ^|(l,ri,0 T ), P 3 = ^(0A-r%O T ), (2) 

respectively, where ta = mA/iriB c and the mass of light pseudoscalar mesons (K, n and 
77 W) has been neglected. For the axial- vector meson A, its longitudinal polarization vector, 
ef, can be defined as 

4 = ^^(1,-^,0 T ); (3) 
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Putting the (light-) quark momenta in B c , A and P mesons as k\, fc 2 , and k 3 , respectively, 
we can choose 



k 1 = (x 1 P+,0,k lT ), k 2 = (x 2 P+, 0,k 2T ), k 3 =(0,x 3 P 3 ~,k 3T ). (4) 

Then, for B c — » AP decays, the integration over fcf , k^, and A;^ will lead to the decay 
amplitudes in the pQCD approach, 



A[B C — > AP) ~ J dxidx2dx 3 bidbib2db2b? J db? J 

•Tr [C(t)<S> Bc (x 1 ,b 1 )<S> A (x2,b2)$p(x 3 ,h)H(x t ,b l ,t)S t (x l )e- s(t) ] (5) 

where bi is the conjugate space coordinate of kiT, and t is the largest energy scale in 
function H(xi, bi, t). The large logarithms \n(mw/t) are included in the Wilson coefficients 



C(t). The large double logarithms (In Xi) are summed by the threshold resummation [23 
and they lead to S t (xi) which smears the end-point singularities on X{. The last term, 
e -s(*) ( i s the Sudakov form factor which suppresses the soft dynamics effectively (24]]. Thus 
it makes the perturbative calculation of the hard part H applicable at intermediate scale, 
i.e., m Bc scale. We will calculate analytically the function H(xi,bi,t) for the considered 
decays at leading order(LO) in a s expansion and give the convoluted amplitudes in next 
section. 

For these considered decays, the related weak effective Hamiltonian H e ff [2l[ is given 

by 

H eff = ^ [V; b V uD (CxMOiCai) + C 2 (ai)O 2 (a0)] , (6) 

with the current- current operators Oi )2 , 

1 = u p -f{l - -y 5 )D a c p >f(l - l5 )b a , 

2 = %7^(1 - 75)^^(1 - lb )b a , (7) 

where V c b, V u d are the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements, "D" denotes 
the light down quark d or s, and Ci(fi) are Wilson coefficients at the renormalization scale 
(jl. For the Wilson coefficients Ci )2 (//), we will also use the leading order (LO) expressions, 



although the next-to-leading order calculations already exist in the literature [21[. This 
is the consistent way to cancel the explicit \x dependence in the theoretical formulae. For 
the renormalization group evolution of the Wilson coefficients from higher scale to lower 
scale, we use the formulae as given in Ref. |5j directly. 



III. ANALYTIC CALCULATIONS IN THE PQCD APPROACH 

In this section, we will calculate the decay amplitudes for 32 charmless hadronic B c — > 
AP/PA decays. Analogous to B c — > PV/VP decays in Ref. jl], there are four kinds 
of annihilation Feynman diagrams contributing to these considered decays, as illustrated 
in Fig. [TJ By analytical evaluation of the two factorizable annihilation (fa) diagrams 
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(a) (b) (o) (d) 

FIG. 1: Typical Feynman diagrams for the charmless hadronic B c — y AP decays. 

Fig. 1(a) and 1(b), we find the corresponding decay amplitude 

F^f = —8TtCpm Bc / dx2<ix3 I b2db2bsdbs 
c Jo Jo 

x {h fa (l - x 3 ,x 2 ,b 3 ,b 2 )E fa (t a ) [x 2 (f>A {x 2 )<f)p(x 3 ) + 2r A r^0p(x 3 ) 

x ((x 2 + 1)^(^2) + (x 2 - 1)^(^2))] + h fa (x 2 , 1 - x 3 , 6 2 , b 3 )E fa (t b ) 

x [{x 3 - 1)0a(x 2 )0£(x 3 ) + 2r A r^ s A (x 2 ) ({x 3 - 2)<j) P P {x 3 ) - x 3 f P {x 3 ))] } , (8) 

where 4> A , 0^' and 4>p' P ' T denote the distribution amplitudes of the axial- vector and 
pseudoscalar mesons, r$ = /m Bc with m p standing for the chiral scale of pseudoscalar 
meson(P), and Cp = 4/3 is a color factor. In Eq. (|S]), the terms proportional to {r A {r p )) 2 
have been neglected because they are small: less than 7% numerically. The function hf a , 
the scales U and Ef a (t) can be found in Appendix B of Ref. [3|. 

For the two nonfactorizable annihilation (na) diagrams Fig. 1(c) and 1(d), all three 
meson wave functions are involved. The integration of 63 can be performed using 5 function 
5(b 3 — 6 2 ), leaving only integration of b\ and b 2 . The corresponding decay amplitude is 

16v^6 f°° 
— 5 — ^Cpm 2 ^ I dx 2 dx 3 / bxdb\b 2 dh 2 

x {Ka( x 2, x 3 ,h, b2)E na (t c ) [(r c - x 3 + 1)4>a(x 2 )4>p(x 3 ) + r A r p (4> s A (x 2 ) 
x ((3r c + x 2 -x 3 + l)(f)p(x 3 ) - (r c - x 2 - x 3 + l)(f>p(x 3 )) + 4> A {x 2 ) 
x ((r c - x 2 - x 3 + l)(f)p(x 3 ) + (r c -x 2 + x 3 - l)(pp(x 3 )))] - E na (t d ) 
x [(r b + r c + x 2 - l)4> A (x 2 )(f)i(x 3 ) + r A r p (0^(x 2 )((4r 6 + r c + x 2 - x 3 
-l)(j)p{x 3 ) - (r c + x 2 + x 3 - l)(j)p(x 3 )) + <f) A (x 2 )((r c + x 2 + x 3 - 1) 
x4>p(x 3 ) - (r c + x 2 -x 3 - l)(f)%(x 3 )))] h d na (x 2 ,x 3 ,h,b 2 )} , (9) 

where r b = mb/rriB c , r c = m c /mB c and r b + r c w 1 for B c meson. 

By exchanging the position of the final state mesons A and P, we can obtain the phe- 
nomenological topology for B c — y PA decays easily. The corresponding decay amplitudes 
for this type of decay channels can be obtained directly by the following replacements in 
Eqs. © and fD), 

<f) A < — y <j)p-, <j) A < — y <f P , <f) A < — > r A < — y r p . (10) 



Before we put the things together to write down the decay amplitudes for the studied 
decay modes, we give a brief discussion about the K\ A -K\p-, fi-fs and hi-h 8 mixing. 



The physical states i^i(1270) and i^i(1400) are the mixtures of the K\a and Kib- 
Kia and K 1B are not mass eigenstates, and can be mixed together due to the strange and 
nonstrange light quark mass difference. The mixing of K\a and Kib can be written as 

1^(1270)) = \K lA )sm0 K + \K 1B )cos6 K , (11) 
|ATi(1400)) = \K 1A )cos8 K - \K 1B )sm6 K . (12) 

If the SU(3) flavor symmetry between (u, d, s) quark was an exact symmetry, K\a and 
Kib would not be mixed with each other. As mentioned in the introduction, the mixing 
angle 9k still not be well determined because of the poor experimental data. In this paper, 
for simplicity, we will adopt two reference values as that used in Ref. jl3j]: 9k = ±45°. 

Analogous to the 77-7/ mixing in the pseudoscalar sector, /i(1285) and /i(1420) (the 
l 3 Pi states) will mix in the form of 

/!(1285)\_/ COS03 sin9 3 \ff 1 \ ( , 

MU20)) \-sm9 3 cos 3 J { fs J 1 J 

Likewise, the /ii(1170) and /ii(1380) (l x Pi states) system can be mixed in terms of the 
pure singlet \hi) and octet \h$), 



(14) 



hi(1170) \ _ f cos 0i sin 0i \ ( h 
/ii(1380) J ~ \- sin 0i cos 0i J \h 8 

where the component of |/i), \hi) and |/ 8 ), \h 8 ) can be written as 

\fi)>\hi) = -4= (I??) + \ss)) , 

\h)Ah 8 ) = ±(\qq)-2\58)), (15) 

where q = (u, d). The values of the mixing angles for l 3 Pi and l x Pi states are chosen as 
13: 



9 3 = 38° or 50°; 9 1 = 10° or 45°. (16) 

By putting all things together, we can write down the general expression of the total 
decay amplitude for the considered decays: 

A(B C -+ AP) = V* b V uD {fsFff^a, + M n r /(PA) d} , (17) 

where a\ = Ci/3 + Ci- Now it is straight forward to present the explicit expressions of 
the decay amplitudes for all 32 considered B c — > AP decays. 
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(1) For AS = processes, 



A{B C ->• 7r + a°) 
„4(B C ->• a+7T°) 



" fB c Fff ai + M^d] } /V2 , 

-^(S c vr+a?) = V^V^ { [fB c F a f f- ai + M^d 

v; b v ud cos { [/ s f;^^ + m:*^] 

+ [/BeF / r i ai + Mr i C 1 ]}/v / 2, 

KtKd sin { [/B C F;r u «i + kt-ci] 
+ [fB C F]r*i+Mr i c 1 ]}/V2, 



.4(£ c ^ 6+tt ) = - 



" [fBF b > ai + M$TC^ } /V2 , 
v4(B c = V; b V ud { [fuFf^a, + M*fici 



v: h v t 

+ 



+ 



^WjJ/^F^ + M^d] 



^(S c ^tt + /i(1285)) 



1 "71" I. fa fa 



sm ' 3 



„4(5 C tt + /i(1420)) 



+F / ?)a 1 + (M:5+M n ^)C' 1 ]} , 

v; b v ud {^[f Bc (F^ + Fjh ai 

) fll + (Mgs + Mft)Ci] } , 



A(B C -+ ^(1170)) = V; h V ud \ ^ [f Bc (F^ + Ff)ai 



V3 

+{M£t+M$)C 1 
-,hU 



+ 



sin^i 



f Bc K 



+Ff) ai + (M^s + M ^) Cl ] } , 
A(B C -+ ^(1380)) = KtK, f/ Bc (< r + 

COS #1 



v/3 



+ 



y/6 



A{B c ^icK l (mQy) = 



AiBc-tlCK^UOO)*) 



V; b V ud {sm9 K [f B Ff a K ^ ai + Mf>-d] 
+ cos^ [f B Ff a K ^ ai + Mf>*d] } , 
V; b V ud [cos9 K [f B FX KlA ai + Mf>^d] 
- sin 9 K [f Bc Ff a K ^ ai + Mj>*d] } , 



+ cosfl* [/a^* V + M^ K d] } , 
^^^(MOO) ^) = V: b V ud [cos0 K [f Bc Fj} AK ai + M^K Cl 

- sin 6 K [f Bc F^ K ai + M^d] } . 

(2) For AS = 1 processes, 

A(B C ^°a+) = ^(S c K+a?) = K* 6 K S {/b c F£° ai ai + Mf>d} 
A(B C ^°&+) = v^(S c = Ok {/i^/fV + M^d} 

.A(fl c ->iiri(1270) 7r + ) = ^(5 c ^i^i(1270)V) 



= O^{sin0* 

r t 

+ cos6 K \ f Bc F fi 



ai 



„4(5 C ->■ iri(1400) 7r + ) = V2A(B C ->• ^(1400)%°) 

r r 7v-n _ 



- sin ^ [/ Bc F£ ? S V + M5 s?r d] } , 
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A(B C -> K + /i(1285)) 



sin 6? 



A(B c ^K + f i(1420)) = 



-2F/f + (Mf/B _ 2M& K )C 1 ] } , 
- sin #3 



^3 



V^6 



„4(5 C ->■ K + /n(1170)) 



,Kh% 
'fa 



-2Fff) ai + (M£* 



2M£*)Ci] } 



a{b c -> K + /n(i380)) = v;;k s 



■ sin 6*i 



+(m^ + m^)c 1 ] + ^[ /Bc (f- 

-2Fff) ai + (M£* - 2M**)C 1 ] } , 



(38) 



*K + (Mf/B _ 2M«*)C 1 ] } , (39) 



(40) 



(41) 



^(5^^(1270)^) = O^{sin0* [/ Bc (cos 0F£ 1Ar " - sin 4>F]f' A \ 

+ (cos0M n x a ^ -sin^M^^d] 
+ COS0* [/ Bc (cos0F^ -sin^F^K 
+(cos^M f £«* - sin^M^^d] } , (42) 

A{B c ^K 1 {lAm) + r 1 ) = V: b V us {cos6 K [/^(cos^F^ - sin^F^K 



+ (cos 
— sin 9 K 



<f>M«^- sin tM^Ci] 



/^(cos^F^-sin^F/f-)^ 



+ (cos0M n x -^ - sin^C^d] } , (43) 
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A(B C -+ K 1 (1270)+7 7 / ) = V; b V us [sm9 K [f Bc (sin cj>Ff a ^« + cos cj>Ff a K ^) ai 

+(sin0M^ + cos^^Cj 
+ cos6 K [/ Bc (sin0F^ + cos cj>Ff a Kw )a x 
+(sin0M*"^ + cos^^Cj } , (44) 

A(B C 1^(1400) V) = V c * b K s {cosfl* [/B c (sin0F^ + <x» <j>Fj' a KlA ) ai 

+(si n 0M^ + cos^^Cj 
-sinfl* [/^(sin^F^ + cos^F^)^ 
+(sin0M^ + cos^lC^d] } • (45) 

IV. NUMERICAL RESULTS AND DISCUSSIONS 

In this section, we will calculate the branching ratios for those considered 32 charmless 
hadronic B c — > AP decay modes. The input parameters and the wave functions to be used 
are given in Appendix El In numerical calculations, central values of input parameters 
will be used implicitly unless otherwise stated. 

For B c — > AP decays, the decay rate can be written as 

T = ^^{l-r\)\A{B c ^AP)f (46) 



where the corresponding decay amplitudes A have been given explicitly in Eqs. 
With the complete decay amplitudes as given in last section, by employing Eq. (T46l) and 
the input parameters and wave functions as given in Appendix |X1 we calculate and present 
the pQCD predictions for the CP-averaged branching ratios of the considered decays with 
errors, as shown in Tables Hl lIVI The dominant errors come from the uncertainties of charm 
quark mass m c = 1.5 ± 0.15 GeV, the combined Gegenbauer moments a* of the relevant 
meson distribution amplitudes, and the chiral enhancement factors tuq = 1.4 ± 0.3 GeV 
and m,Q = 1.6 ± 0.1 GeV, respectively. 

Based on the numerical results as given in Tables I -IV, we have the following remarks: 

• The pQCD predictions for the CP-averaged branching ratios of considered B c decays 
vary in the range of 10~ 6 to 10~ 8 . There is no CP violation for all these decays within 
the standard model, since there is only one kind of tree operator involved in the 
decay amplitude of all considered B c decays, which can be seen from Eq. (fTj 



Among the considered B c — > AP decays, the pQCD predictions for the branching 
ratios of those AS = processes are generally much larger than those of AS = 1 
channels (one of the two final state mesons is a strange meson), the main reason is 
the enhancement of the large CKM factor (Vud/Vusl 2 ~ 19 for those AS = decays 
as generally expected. For B c — > (af, bf) (tt°, K°) decays, however, the difference is 
not so large, because the enhancement due to the CKM factor is partially canceled 

by the differences between the magnitude of individual decay amplitude \Fj a 1 
and|F;5*°|. 
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For B c — > (ai,bi)n decays, the same component of uu — dd is involved in both 
axial- vector (a®, 6°) and the pseudoscalar ir° meson at the quark level. We therefore 
find the same branching ratios for B c — > n + a\ and B c — > a^7r°, and for B c — > n + b\ 
and B c — > bfir , respectively. 

From the numerical results as shown in Table HI one can see that: 

Br(B c 6itt) ~ 14 x Br(B c ai ir), 

Br(B c -> b x K) ~ 16 x Sr(S c ->■ Oilf). (47) 



This pattern agrees well with that as given in Ref. [14j . Il6 |. 
Unlike B c — >■ (ai, &i)(vr, iT) decays, we find that 

Br(B c -> a 1 (7 7 ,7 7 / )) ~ 5r(£? c -> ( 48 ) 

The main reason is that the suppressed factorizable annihilation amplitudes cancel 
each other for B c — > a\T]^' decays, while the enhanced nonfactorizable ones cancel 
each other for B c — > bir)^ decays. 

For B c ->■ Z°(ir 1 (1270) + ,/^i(1400) + ) and B c -»■ (#i(1270)°, Ki(Umf)K+ decays, 
their BRs strongly depend on the value of the mixing angle Ok of the Kia-K 1b 
system. From Table II, one can see that 

< 2, (49) 



Br(B c iTKi(1400) + ) Br(B c ->• if +if i(1400) c 



£r(B c -> lT#i(1270)+) ^(5 C K+Ki(1270)°) 
for = 45°, while 

Br(B c -»■ 7rVi(1400) + ) _ 5r(5 c -»■ K+i?i(1400)°) _ 1 
fir( J B c ^K°A^ 1 (1270)+) ~ Br(B c -» ^+^(1270)°) ~ 2' (5 ° } 

for 0^- = —45°. This means that one can determine the sign and size of 9 k after 
enough B c events become available at the LHC experiment. 



TABLE I: The pQCD predictions of branching ratios(BRs) for B c — > (a\,b\)P decays. The 
source of the dominant errors is explained in the text. 



AS = 
Decay modes 


BRs(l(T 7 ) 


AS = 
Decay modes 


BRs(10" 6 ) 


B c -»• 7T+a? 
B c -»• a^vr 


S-OigiW^^aOiiKmo) 
2.9«i(m c )+?; 2 7 (a,)l^(m ) 


B c vr+b? 
B c b+^° 


4.3llj(m c )li|(a i )l°;?(m ) 


B c ->■ &fr] 
B c -> a+r/' 


6.8^ : 4 (mc )+|r (a4 )+o : o (mo) 
4.6i^(m c )t};I(a i )iH(^o) 


B c b+7? 
B c b+7?' 


0.6lg5(m 8 )lg5(a i )lH(mo) 
0.4i^(m c )^i(a,)T°;°(m ) 


AS = 1 
Decay modes 


BRs(l(T 8 ) 


AS = 1 
Decay modes 


BRs(l(T 7 ) 


B c -> a+K° 
B c -»• K+a? 


3.4l^(m c )iO( aj) +o.6 (mo) 

i-7i^K)li:?(ai)ig:?(m ) 


B c b+K° 
B c K+b° 


5.4t°|(m c )+|2 (aj) +o.2 (mo) 
2.7^^(m c )l^( ai )^(mo) 
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TABLE II: Same as TableHbut for B c (Ki(1270),Kx(U00)){Tr,K,r),rj') decays. 



AS = 
Decay modes 


BRs(10~ 7 ) 
^ = 45° 


BRs(lCT 7 ) 
9 K = -45° 


B c -)• K U Ki(1270)+ 
B c -> K K!(1400)+ 


8 2 + Hfra r ) + J 6 , 3 fe) + S'5(mn) 
17.3l^(m c );?^(a,);?;°(m ) 


1 •^ t — 4.1 \""cj — 16.1 V"*/ — 1.5 V" Hjy 


B c _> K 1 (1270)°K + 
B c ^ Ki(1400)°K+ 


15.8 + ri(m r ) + J B ; 6 fa,-) + nSfmn) 

-3.31 w-8.1 V 1 ^^/— O.OV U/ 

31.7^ 2 3 (m c )l?° : °(a,)lS;«(mo) 


32.0 + ^ 4 (m c ) + ?Q7(a i ) + °-°(m ) 

" 7.3 V c/— 19.7V * /—-2.0V u / 


AS = 1 
Decay modes 


BRsflO -8 ! 
= 45° 


BRsflO -8 "! 

e K = -45° 


B„ ->■ Ki fl270)°7r + 
B c ^Ki(1400) vr+ 


U, °-3.3V c /— 4.5 V"*/ — 1.3 V'-U; 

5.8i^(m c )i^( ai )iH(m ) 


xj - v — o.7\ ln cJ — i.gvu-ij-o.l 

6.8^;°(m c )l^( aj )+o : 7 (mo ) 


B c ^Ki(1270)+^° 
B c ^Ki(1400)+7T 


3.4 + ^( mc )+3-3( a .)+0-4( mo ) 

2.9l°i(m c )ljj(a i )^(m ) 


3.0 + °- 7 (m c )+}' 7 (a 4 ) + [Ji(mo) 
3.4^ : 5 7 (m c )i^(a,)i° : 4 7 (m ) 


AS= 1 
Decay modes 


BRs(10" 8 ) 
^ = 45° 


BRs(10" 8 ) 
^ = -45° 


B c ->• Ki(1270)+77 
B c -»• Ki(1400)+77 


16.8^;°(m c )t^.i (aj) +o.o (mo) 
26.9i|l(m c )Ht 8 7 (a i )lH(^o) 


27.2^ : °(m c )+^; 8 (a l )+°;°(mo) 


B c -»• Ki(1270)+7?' 
B c ^Ki(1400) + 7?' 


2.7i^(m c )i|2 (ai )+o.o (mo) 


n.6lil(mc)lS5(«i)i8:8(mo) 

2-7l^(m c )l^(a i )lH(m ) 



TABLE III: Same as TableHbut for B c -> (/i (1285), fx (1420)) (tt, if) decays. 



AS 1 = 
Decay modes 


BRs(10" 8 ) 
fl 3 = 38° 


BRs(10" 8 ) 
3 = 50° 


B c -> 7r+fi(1285) 
B c -»• 7r+fi(1420) 


52.4i 9 7 i( mc )i30 : 2 (aj )+2i.7 (mo) 
8.5;?i(m c )^ : o (ai) +o.r (mo) 


44.8l 7 ;° 7 (m c )^.6 K)+ i9.8 (mo) 
16.0i| 8 (m c )« 1 2 5 (a i )l?|(m ) 


AS = 1 

Decay modes 


BRs(10" 8 ) 
fl 3 = 38° 


BRs(10" 8 ) 
3 = 50° 


B c K+fi(1285) 
B c -> K+fi(1420) 


1.6l ;° 7 (m c )t^(a,)^i(m ) 
7-4i^(m c )t||(a i )i^(m ) 


1.5i^(m c )^(a,)i^(m ) 
7.5t^(m c )t|2 (ai) +0.4 (mo) 



For the AS* = 1 B c — > Kyn decays, their decay rates have a very weak depend 
on the value of mixing angle Qk'- 

Br(B c ir 1 (1270) Tr + ) « Sr(5 c -+ #i(1400)°7r + ) « 6 x 1(T 8 , 
£r(£ c -> K 1 (1270) + 7T°) « Br(B c -> Ki(1400) + tt°) « 3 x 1(T 8 , 

for both 6^ = 45° and —45°. This point will also be tested at LHC. 
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TABLE IV: Same as TableEQbut for B c -»• (hi (1170), hi (1380)) (tt, X) decays. 



AS = 
Decay modes 


BRs(10" 8 ) 
01 = 10° 


BRs(10" 8 ) 
01 = 45° 


B c -> 7r+hi(1170) 
B c -»■ 7r+hi(1380) 


60.5i^ (r n c )t32.i (aj)+ 25 o 8 (mo) 
1.2^ : ^(m c )t^ (ai) +o.5 (mo) 


49.llf 9 8 (m c )^; 9 7 (a i )^ 5 7 (m ) 
12.4i3 : 9( mc )+3.3 (aj) +5.5 (mo) 


A5 = l 
Decay modes 


BRs(10" 8 ) 
#i = 10° 


BRs(10" 8 ) 
01 = 45° 


B c -»• K+hi(1170) 
B c -> K+hi(1380) 


14.9l?;°(m c )^ 2 6 (a i )^(m ) 
22.2i^ 3 5 (m c )li^(a i )l°;°(mo) 


16.8l^(m c )^ : 9 (ai) +o.3 (mo) 
20.2^° 5 6 (m c )^ 1 2 7 (a l )tH(^o) 



the pQCD predictions have a strong Q K dependence: 



• For B c Ki7]® decays, 



Br(B c -> ifi(1400) + r?) 
£r(5 c ->• ATi(1270)+7/) 
Sr(S c -> iTi(1400) + ^) 



Br(B c -> iri(1400) + r/) 1 



(53) 



Br(B c -> Ki(1270)+V) 

for 0# = 45°, while 



Sr(S c -> ^(1270)+^) 1.6' 
gr(B c -»^i(1400)V) ^ J_ 

Br(B c ^ Kxifflfy+Tf) ~ 4.3' 1 J 

for #k = —45°. It is easy to see that these B c — > Kirj^ decays are sensitive to the 
mixing angle 9k- Analogous to the B c — > K*(rj, if) decays the above four decays 
are dominated by the factorizable annihilation diagrams. 

The theoretical predictions for the branching ratios of B c — > Ki(1270)P and B c — > 
Ki(1400)P for 9k = 45°, as listed in the column two of Table II, are roughly 
exchanged with respect to those of the third column for the choice of 9k = —45°. 
Such simple relation comes from the fact that the two states 7£i(1270) and _Kd(1400) 
can go one into another as a mixture of K X a and Kib states when one sets the mixing 
angle 9k = 45° or —45° respectively, as can be seen from Eqs. (11,12), 

|#l(1270)> | 0K= 45° = |^l(1400)>| %= _45°, 

|#i(1400)>| -|Xi(1270))| (55) 

This relation further leads to the following relations between the decay amplitudes 
of B C ^K X P: 

A(B c ^K 1 (1270)P)\ 6k = 45 o = A(B C ^ K 1 (U00)P)\ dK= _ 45 o , 

A(B C -> Ki(1400)P)|, x=45 o = -A(B C -»• K 1 (1270)P)\ 9k =-4 5 ° , (56) 
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and finally we obtain the special pattern of branching ratios as listed in Table II. 
The small differences in corresponding decay rates are due to the difference in the 
masses of lfi(1270) and i^i(1400) mesons. The numerical relations as shown in 
Eqs. (I49ti50p and (I53ti54|) are also induced by the same mechanism. 



For the four B c — > fx(K, tt) decays, one can see from Table HTT1 that 
Br(B c -> tt+/i(1285)) _ / 6.2 for 6 3 = 38° 

and 



Br(B c -> tt+/i(1420)) i 2.8 for 6 3 = 50° 
Br(B c -> ir+A(1285)) 



(57) 



Br(B c -> AT+/i(1420)) 



0.2 (5* 



for # 3 = 38° and 50°. The relations in Eqs. ( I57|5"8|) can be understood as follows: 
(a) since /i(1285) and /i(1420) are the mixed states of fi and /g (see Eq. f JTB"]) ) and 
both sin #3 and cos 6*3 are positive for #3 = 38° and 50°, the contribution from the 
common component (qq) of fi and f$ will interfere constructively (destructively) 
for B c — > 7r + /i(1285) (B c — > 7r + /i(1420)) decay, this results in the large difference 
for the decay rate of the two decays; (b) for the two AS = 1 decays, however, 
the new component (ss) will provide additional contributions to the considered 
decays. Furthermore, the contributions from (ss) and qq interfere constructively for 
B c -> K+f x (U20), but destructively for B c -> K+ ^(1285) decay. 

• The pQCD predictions for B c — > (/i 1 (1170), hi(1380)(K, tt) decays, as given in Table 
IV, can be explained in a similar way as for B c — >■ (/x(1285), /i(1400)(X, tt) decays. 

• Since the LHC experiment can measure the B c decays with a branching ra- 
tio at 10 -6 level [2|, our pQCD predictions for the branching ratios of B c — > 
K(Ki(l270),Ki(l400)) and biix decays could be tested in the forthcoming LHC 
experiments. 

It is worth stressing that the theoretical predictions in the pQCD approach still have 
large theoretical errors induced by the still large uncertainties of many input parameters, 
e.g. Gegenbauer moments at. For most considered pure annihilation B c decays, it is hard 
to observe them even in LHC due to their tiny decay rate. Their observation at LHC, 
however, would mean a large non-perturbative contribution or a signal for new physics 
beyond the SM. 

We here calculated the branching ratios of the pure annihilation B c — > AP decays 
by employing the pQCD approach. We do not consider the possible long-distance (LD) 
contributions, such as the re-scattering effects, although they may be large and affect the 
theoretical predictions. It is beyond the scope of this work. 

V. SUMMARY 

In short, we studied the charmless hadronic B c — > AP decays by employing the pQCD 
factorization approach based on the kr factorization theorem. These considered decay 
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channels can occur only via the annihilation diagram in the SM and they will provide 
an important platform for testing the magnitude of the annihilation contribution and 
understanding the content of the axial- vector mesons. 

The pQCD predictions for CP-averaged branching ratios are displayed in Tables (TTl fTVT) . 
From our numerical evaluations and phenomenological analysis, we found the following 
results: 

• The pQCD predictions for the branching ratios vary in the range of 10~ 6 to 1CT 8 . 
The B c ->■ K°(K 1 (1270) + ,i ; s:i(1400) + ) and other decays with a decay rate at 1(T 6 
or larger could be measured at the LHC experiment. 

• For B c — > A P decays, the branching ratios of AS = processes are generally much 
larger than those of AS = 1 ones. Such differences are mainly induced by the CKM 
factors involved: V u d ~ 1 for the former decays while V us ~ 0.22 for the latter ones. 

• Since the behavior for x Pi meson is much different from that for 3 Pi meson, the 
branching ratios of pure annihilation B c — > y4( 1 P 1 )P are basically larger than that 
of B c A( 3 Pi)P, which can be tested in the LHC and Super-B experiments. 

• The pQCD predictions about the branching ratios of B c — > Kii]^ and K\K decays 
are rather sensitive to the value of the mixing angle 9 k- One can determine 9 k 
through the measurement of these decays if enough B c events become available at 
the LHC experiment. 

• The pQCD predictions still have large theoretical uncertainties, mainly induced by 
the uncertainties of the Gegenbauer moments a« in the meson distribution ampli- 
tudes. 

• Because only tree operators are involved, the CP-violating asymmetries for these 
considered B c decays are absent naturally. 
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Appendix A: Input parameters and distribution amplitudes 



The masses (GeV), decay constants (GeV), QCD scale (GeV) and B c meson life- 
time (ps) are 



A (^4) 
MS 


= 0.250, 


m w 


= 80.41, 


m Bc = 


6.286, 


/s c = 


0.489, 


m ai 


= 1.23, 


fax 


= 0.238, 


m bl = 


= 1.21, 


fax - 


0.180, 


m K 1A 


= 1.32, 


fK 1A z 


= 0.250 


m K 1B = 


= 1.34, 


fK 1B 


= 0.190 


ffx 


= 0.245, 


m h 


= 1.28, 


fh = 


0.239, 


m h = 


= 1.28, 


fax 


= 0.180, 


m hl 


= 1.23, 


fas = 


0.190, 


m hs -- 


= 1.37, 


ml 


= 1-4, 


< 


= 1.6, 


Vq 

TUq = 


1.08, 


m v s = 


1.92, 


m b 


= 4.8, 


fa 


= 0.131, 


f K = 


0.16, 


T Bt = 


0.46 . 



(Al) 

For the CKM matrix elements, here we adopt the Wolfenstein parametrization for the 
CKM matrix, and take A = 0.814 and A = 0.2257, p = 0.135 and fj = 0.349 pi . 

For the distribution amplitudes of pseudoscalar mesons, we adopt the same forms as 
that used in the literature (See Ref. |3j and references therein). 

The twist-2 distribution amplitudes for the longitudinally polarized axial-vector 3 P X 



and x Pi mesons can be parameterized as [13l . Il8 

/ 



(p A (x) 



6x(l 



x 



2^/2N ( 

As for twist-3 LCDAs, we use the following form: 



Oq + 3a| t + a\ ~(5t 2 



1) 



<f>\{x) 



3/ 



2^2N~ C 
3/ 



d 



a^t 2 + -ait(3t 2 - I) 
{x(l - x)(a£ + a±t)} 



(A2) 

(A3) 
(A4) 



2^/2N c dx 

where / is the decay constant and t = 2x — 1. It should be noted that for the distribution 
amplitudes of strange axial- vector mesons and Kib, x stands for the momentum 
fraction carrying by s quark. 

Here, the definition of these distribution amplitudes (ft a (x) satisfy the following relation: 



iPl [X) 



II 3 fi 

where we have used all' =1. 



2v^AT 



i Pl {X) 



2^2W C ' 



(A5) 



The Gegenbauer moments have been studied extensively in the literatures (see Ref. [13 
and references therein), here we adopt the following values: 



11,11 


= -0.02 ±0.02; 


_L,ai 


= -1.04 ±0.34; 


af 1 = 


-1.95 ±0.35; 




= -0.04 ±0.03; 




= -1.06 ±0.36; 


1 l,>u 

a{ = 


-2.00 ±0.35; 




= -0.07 ±0.04; 




= -1.11 ±0.31; 


4 M = 


-1.95 ±0.35; 




= 0.00 ±0.26; 




= -0.05 ±0.03; 


-L,Kxa 


0.08 ±0.09; 




= -1.08 ±0.48; 


a ll ' KlB 


= 0.14 ±0.15; 


ap lfl = 


-1.95 ±0.45; 




= 0.02 ±0.10; 




= 0.17 ±0.22. 







(A6) 
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